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Abstract. - The formation of normal-state domains in type-I superconducting indium films 
is investigated using the high resolution magneto-optical imaging technique. The observed 
patterns consist of coexisting circular and lamellar normal-phase domains surrounded by the 
superconducting phase. The distribution of domain surface areas is found to exhibit a threshold, 
above which only the lamellar shape is observed. We show that this threshold coincides with the 
predicted critical surface area for the elongation instability of the circular shape. The partition 
of the normal phase into circular and lamellar domains is determined by the combined effects 
of the elongation instability and the penetration of magnetic flux by bursts at the early stage 
of pattern formation. It is not governed by mutual interactions between domains, as usually 
assumed for self-organized systems. 



Introduction. - A spontaneous phase separation into domains is encountered in a large 
number of systems including magnetic fluids [1-3] , Langmuir monolayers confined at air- water 
interface [4], ferro- and fcrrimagnetic layers [5,6], adsorbates on a metal substrate [7], and 
type-I superconductors in the intermediate state (IS) [8]. The formation of domains originates 
from the balance between the short-range attractive interaction associated with the interfacial 
tension between the two phases and long-range interactions. The regular shapes of the domains 
are generally not stable due to the effect of long-range interactions: lamellar domains present 
undulation or peristaltic instabilities [3,5,9-11], circular domains (bubbles) may elongate and 
produce finger structures [2,4,12,13]. A crucial point for understanding the dynamics of 
domain patterns is to determine the contribution of these instabilities to the morphogenesis 
of domains [3,14-18]. 

In type-I superconductors, domain patterns are observed in film samples submitted to 
a perpendicular magnetic field. The IS pattern consists of coexisting normal-state (NS), 
flux-bearing domains, and superconducting (SC) domains. By analogy with magnetic fluids, 
the instability of bubble domains was proposed as the mechanism of formation of ramified 
structures [15]. However, the progressive transformation of a circular domain into a ramified 
structure, which is encountered in ferromagnetic films [12] and in magnetic fluids [2], has 
not been reported in type-I superconductors up to now. This raises the question of how 
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unstable NS bubbles can be produced in SC systems. The stability limit for the circular 
bubble depends on its diameter and on the magnetic Bond number Bm that characterizes the 
ratio of the magnetic energy and the interface energy [14]. Tuning the external field changes 
Bm in magnetic fluids [14] and the diameter of bubbles in ferromagnetic systems [6]. This 
opens the possibility to drive the system sufficiently out of equilibrium to bring bubbles above 
their stability limit. In contrast, for superconductors, increasing the external field leaves 
unchanged both Bm and the diameter of the NS bubbles in the IS pattern [19,20]. 

Bubble instability could occur during the penetration of the magnetic flux into the su- 
perconductor. In type-I superconductors, the early penetration of magnetic flux proceeds by 
magnetic flux bursts from the sample edges [21-24] . Studies of bursts remain beyond the scope 
of experimental investigation due to their high propagation velocity (w Ini/s) [22]. There- 
fore, little is known about their characteristics. Whether the domains produced by the bursts 
are sufficiently out of equilibrium to overstep the bubble stability limit remains a fully open 
question. 

This letter presents an analysis of the volume and shape of NS domains. A clear correlation 
between the surface area and the shape of domains is revealed. The threshold area, above 
which no more cylindrical domains are observed, is shown to correspond to the elongation 
instability limit. We discuss the connection between this instability and the penetration by 
flux bursts. Finally, the contributions of these processes to the formation of domains in type-I 
superconductors are emphasized. 

Observation of domains in type-I superconductors. - The IS flux pattern was studied 
by high-resolution Faraday microscopy imaging, a technique that probes the normal compo- 
nent of the local induction at the top surface of a superconductor. Experimental details are 
given elsewhere [20]. Indium fllms with thickness d =2.2 and 10.0 ± 0.1 /im were grown by 
evaporation directly onto a magneto-optic layer [25]. A 112 ± 1 fim thick In sample was cut 
out from a GoodFellow 99.9 % purity annealed foil and placed against a magneto-optic layer. 
The critical temperature and critical field for indium are Tc^S. 41 K and Hc{T—0)—28.2 mT, 
respectively. Hc{T) was assumed to follow a Bardeen-Cooper-Schrieffer temperature varia- 
tion: Hc{T) = i7c(0)(l — (T^/T^)). The samples were immersed into superfluid helium at 
temperature T < 2 K. They were first zero-field cooled then subjected to field cycles. The 
external field Hq was applied perpendicularly to the sample. 

Figure n shows two typical IS patterns obtained close to the edge of the sample when the 
field is increased. When the reduced apphed field h = Ho/Hc{T) is sufficiently large, the 
density of domains is homogeneous even close to the edge (fig. ft,=0.345). This is not 
the case at low field value (fig. h— 0.105) where the domain pattern is separated from 
the sample edge by a full diamagnetic band about 50 //m wide. This diamagnetic band is 
the signature of the geometrical barrier that controls the penetration of magnetic flux. This 
penetration from the edge to the inner part of the sample occurs sequentially, by bursts, each 
carrying a small amount of flux [21-24]. 

The two typical shapes of NS domain are clearly seen in flg. ^ almost circular domains 
(bubbles) and lamellar domains. The shape of domains is usually assumed to be determined 
by the interface energy and the interactions between domains, as for example in ferromagnetic 
films [6]. However, in superconductors, the penetration of magnetic flux by bursts may play 
a predominant role in the first stage of formation of IS patterns. In order to get more insight 
into this point, we performed a statistical analysis of the NS domain shape. 

Analysis of domain shapes. ~ As can be qualitatively observed in fig.Q circular domains 
seem to have a smaller surface area than lamellar ones. In order to discriminate between 
bubbles and lamellae in a quantitative manner, the domain perimeter p was systematically 
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Fig. 1 - Typical IS patterns in an indium sample with thickness d—10 /im for increasing values of the 
applied magnetic field. Images (a) and (b) correspond to ft=0.105 and h— 0.345, respectively (T =1.85 
K). Flux-bearing NS domains appear in black and have circular or lamellar shapes. The edge of the 
sample is along the right edge of the image. The few white domains correspond to magnetic flux 
which was trapped at zero field [25]. 



measured as a function of their surface area A. Typical results obtained for different values 
of h are presented in fig. El for the 2.2 /im thick sample. The maximum value of A increases 
noticeably with h. The weak dispersion of data indicates that p and A are strongly correlated. 
Moreover, the data points, shown in a log-log plot, are essentially grouped along two straight 
lines. For the lowest values oi A (A < 100 fJ.m^), the solid lines in fig. |2] correspond to the 
equation of a circle p — 2%/ nA. The agreement with the experimental data simply indicates 
that domains with small area have a circular shape. The dashed lines plotted in fig. [3 for 
large areas (A > 100 ^m^) correspond to the equation p = 2A/w, which is valid for long 
lamellae. For these plots, w is the lamella equilibrium width calculated in the framework of 
the "constrained current- loop" model [20]. The good agreement with the experimental data 
shows that large area domains have a lamellar shape. We conclude that the shape of a domain, 
either circular or lamellar, is determined by its surface area. 

Moreover, as indicated by arrows in fig. |21 there is a threshold value for A, above which 
the circular shape is absent from the IS pattern. We define precisely this threshold Ath as the 
value of A above which no more data point is found on the line p = \/4ttA. We define here 
the reduced threshold diameter of bubbles 2Rth/d, where Rth = \/ Athji^- The variation of 
IRthjd with h is shown in fig.^lfor sample thicknesses d =2.2 and 10.0 [ivt\. When h is raised 
from to 0.4, the volume fraction of the sample occupied by NS domains increases from to 
0.5. However 2Rth/d is found to increase by only 20% and 10% for the 2.2 and the 10.0 /im 
samples, respectively. This means that 2Rth/d only weakly depends on mutual interactions 
between domains, and suggests a comparison of the threshold diameter with the stability limit 
of an isolated circular domain [2, 14]. 

Stability limit for the circular shape. - In order to determine the critical diameter above 
which the circular shape becomes unstable, we consider the linear dynamics near the circular 
initial condition. The calculation is performed in the framework of the "constrained current- 
loop" model [20]. Let us assume that the applied magnetic field Hq produces a set of N NS 
circular domains of radius r, each bearing an induction /i„. The constraint of flux conservation 
is written Hq = p„/i„ where p„ = Nnr^d/V is the volume fraction of the normal phase in 
a sample of volume V. In the limit of non-interacting bubbles {Hq —>0), the free energy per 
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Fig. 2 - Typical correlation between the perimeter p of the domains and their surface area A, shown 
here for a 2.2 /j,m thick indium sample (T =1.815 K). The data, reported for three values of the 
increasing applied magnetic field h = Ho/ He are displayed in a log-log plot. For A <1Q0 /xm'^, the 
solid lines correspond to p = 2\/^jvA, i.e. to circular domains. For A >100 /xm^, the dashed lines 
correspond to p = 2A/w, i.e. to lamellar domains much longer than wide, where w is the calculated 
lamella width {Bm =0.73). The arrows indicate the threshold area value for which the domain 
shape starts to deviate from the circular one. The critical area for the elongation instability A"^^ is 
indicated by the dotted line. 



unit volume of the system is 
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where Fg is the free energy density of the SC state. The two terms in the first parenthesis 
are the "bulk" terms corresponding to the condensation energy and the magnetic energy 
respectively. The first term in the second parenthesis is the interface energy where the surface 
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Fig. 3 - Threshold diameter 2Rth/d (black squares) as a function of the applied magnetic field 
h = Hq/Hc obtained for the 2.2 /im (a) and the 10 ^m (b) thick indium samples. The empty circle, 
with the dotted line as a guide for the eye, corresponds to the critical diameter 2Rc/d for isolated 
circular domains predicted by eq. |1] {Bm = ci/27rA(T)=0.73 for d=2.2 fim and 3.2 for ci=10 ^im, 
respectively. A(T) is assumed to follow the empirical law A(T) = A(0)/^1 - (T/Tc) with A(0) = 
0.33 fj,m [26]). The dashed lines correspond to the equilibrium diameter of circular domains calculated 
at the transition from the bubble hexagonal phase to the lamellar phase. The line length represent 
the field range in which the energies of the two structures differ by less than 0.3 %. 
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tension apfs is equal to AH'^/Stt, with A the interface energy parameter. The last term is 
the correction of the magnetic energy due to the curvature of the flux lines in the free space 
on both sides of the sample. It can be calculated using the Biot-Savart interaction of the 
screening current encircling NS domains [15,20]. Minimizing F with respect to p„ and r 
yields hn = He and the equilibrium bubble diameter. 

For the calculation of the stability limit, the domain wall motion is supposed to result from 
the balance between the generalized force er6F/6r and a local viscous force —rjerdr /dt [14,15]. 
The resulting equation of motion is 

-^e.| = .a^s + ^rer, (2) 

where k is the interface curvature. Consider small deviations ^ in the radius of the circle, 
parametrized by the polar angle ip: r{ip, t) — [Rq +£,{ip,t)] er{ip) and suppose that the per- 
turbation is in the form ^{ip, t) « exp{ant)cos{nip) . The growth rate of mode n is then given 

by 

an^[{l~n^)/Rl + 8B„,/3d^]aNs/v , (3) 

where the magnetic Bond number i?,„ is equal to (i/27rA(T). Note that the same result is 
obtained for magnetic fluids in the limit d oo [3, 14]. The stability limit is reached for 
CTn = 0: 

2Rl^ / 3(n2-l) 
d Y 2B„ ^ ' 

The reduced critical diameter 2i?"=^/c? predicted for the elongation instability is plotted in 
fig- Eland compared to the threshold diameter 2Rth/d, measured for the 2.2 and the 10.0 /im 
thick samples. Within the experimental error bars, 2i?"^^/d coincides with the extrapolation 
of 2Rth/d for ft, — > 0. The variation of the stability limit with increasing inter-domain interac- 
tions remains beyond the scope of this paper [27]. However, one can compare the maximum 
threshold, 2Rth/d at /i w 0.4, to the predicted equilibrium bubble diameter for which the free 
energies of an hexagonal array of bubbles and a ID-array of lamellae are equal (fig. 01 • The 
good quantitative agreement indicates that the measured variations of the threshold 2Rth/d 
with h are consistent with the two limits for the critical diameter: for h close to zero or close 
to the bubble-lamellae transition field. Moreover, 2Rth/d(h 0) was systematically mea- 
sured for different sample thicknesses. The results are shown in fig. 0] The same figure also 
displays the result obtained for Hg in ref. [28] for p„ « 0.3 - 0.4 with A(0) = 0.084 /im [19]. 
All the data gather on a single master curve indicating that 2Rc/d and _B,„ are appropriate 
reduced variables. A very good quantitative agreement is obtained with the prediction of eq. 
01 without any adjustable parameter. Note that the magnetic Bond number Bm is varied by 
more than two orders of magnitude (0.72-150). This demonstrates that the shape (circular 
or lamellar) of the NS domains is controlled by the elongation instability limit of the circular 
shape. 

Discussion. - This result raises the general question of the role of the elongation insta- 
bility in the formation of patterns in type-I superconductors. It was proposed in Ref. [15] 
that the formation of lamellar structures could originate from the instability of bubbles. A 
careful examination of NS domains in the 2.2 and the 10.0 /im thick samples indicates that 
the NS bubbles start to disappear only when h reaches ~ 0.3 — 0.4. Above this limit, the dis- 
appearance of bubbles is found to result from their fusion with other bubbles or with already 
existing lamellae. Below /i « 0.3 — 0.4, the increase of h does not change NS bubble surface 
area. This is a consequence of the conservation of magnetic flux ($ « HcA) in a NS domain 
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Fig. 4 - Reduced critical bubble diameter 2Rc/d as a function of the magnetic Bond number Bm 
shown in a log-log plot. The filled circles represent the threshold diameter in the limit h — »0 for In 
samples. The empty circle is reported from ref. [28] for Hg. The shaded region corresponds to type-II 
superconductivity for very thin samples. The solid line represents the critical diameter of an isolated 
bubble predicted by eq. |1] A good agreement is obtained over more than two orders of magnitude of 
Bm without any adjustable parameter. 

isolated in the SC matrix [19]. Therefore A cannot increase beyond the stability limit. We 
conclude that, unlike magnetic fluids [2] and ferromagnetic films [12], the early formation of 
lamellar structures in type-1 superconductors does not result from the instability of bubbles 
already present in the IS pattern. 

We now address the question of the contribution of flux bursts to the shape of NS domains. 
Isolated NS domains appearing at low fleld are produced by magnetic flux bursts. The mech- 
anisms of production, as weU as the distribution of amplitude of the bursts, remain to a large 
extent unexplored. When h is increased from zero, NS domains first appear as bubbles. With 
further increase of the field, lamellar domains appear [20]. As seen in fig.|3 this leads, for the 
2.2 /im sample, to a large increase of the fraction of lamellae between h — 0.067 and h = 0.2. 
Those lamellae do not arise from the fusion of bubbles since the latter start to disappear only 
when h reaches « 0.3 — 0.4: they are directly produced by flux bursts. The penetration of 
magnetic flux bursts thus results in the production of both bubbles and lamellae. 

Furthermore, images of the IS pattern show a leveling-off of the number of NS bubbles 
produced by flux bursts when the external field is increased [20]. For the 2.2 and the 10.0 /im 
thick films, this leveling-off occurs at /i = 0.10 — 0.15 while the free energy of the bubble lattice 
remains lower than the free energy of the lamella lattice up to /i « 0.3 [20]. As a consequence, 
the leveling-off of the number of NS bubbles is not governed by the competition between the 
interface energy and long-range magnetic interactions, as assumed in equilibrium mean-field 
models [18,20]. 

In view of the above experimental results, we propose the following scenario for the forma- 
tion of NS domains. The shape of NS domains is determined during their penetration through 
the diamagnetic band. The amplitude of flux bursts, i.e. the surface area A = ^/ He oi the 
NS domains released on the edge of the sample, increases with h. For small values of /i, the 
released domains have a surface area smaller than the critical area for the elongation insta- 
bility A"^^: they become bubbles. With further increase of h, the surface area of the largest 
released domains exceeds they become lamellae. The appearance of new bubbles stops 

for a field value h such that the surface area of all the released domains is larger than . 
Therefore the NS domain shape depends on the distribution of amplitudes of flux bursts and 
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its evolution with the external field with respect to the critical area for the elongation insta- 
bility. We emphasize that the process proposed here for domain formation concerns only the 
low field range. With increasing field, NS domains become connected to the sample edge and 
the growth of the normal phase proceeds by continuous penetration of the magnetic flux [20] . 

Conclusion. In conclusion, our study shows clear evidence of the existence of the elon- 
gation instability of the circular shape in type-I superconductors. The role of this instability 
in the formation of domain pattern is clarified. At the early stage of fiux penetration, the par- 
tition between bubbles and lamellae is not driven by a thermodynamical equilibrium resulting 
from the competition between long-range and short-range interactions. It is determined by 
both the bubble elongation instability limit and the distribution of amplitudes of fiux bursts 
penetrating on the edges of the samples. These results explain well why almost perfect hexag- 
onal arrays of bubbles are hardly observed in type-I superconductors while they are commonly 
reported in ferromagnetic films and magnetic fiuids. We hope that this work will stimulate 
further theoretical work in order to get more insight into the dynamics of penetration of the 
magnetic flux. 

* * * 

The authors gratefully acknowledge stimulating and fruitful discussions with A. Cebers. 
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